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Abstract.  This  paper  presents  numerical  results  obtained  by  applying  an  LES  model  based  on  a 
rational  (subdiagonal  Pade)  approximation  in  the  wave  number  space.  We  used  a spectral  element 
code  to  test  this  LES  model  and  the  Smagorinsky  model  with  Van  Driest  damping  in  the  numerical 
simulation  of  a 3D  channel  flow  at  ReT  = 180.  The  corresponding  results  were  compared  with  the 
fine  DNS  simulation  of  Moser,  Kim,  and  Mansour. 


1.  Introduction 

Introduced  in  1970  by  Deardorff  [4],  Large  Eddy  Simulation  (LES)  is  one  of  the  most  promising 
approaches  to  the  numerical  simulation  of  turbulent,  flows.  LES  has  a rich  history,  being  success- 
fully developed  and  applied  to  a wide  range  of  applications  by  the  engineering  and  geophysics 
communities. 

LES  is  based  on  the  idea  that  the  large  scales  in  the  flow  are  essential,  being  responsible  for 
the  important  engineering  quantities  (heat  transfer,  mixing,  etc.),  whereas  the  small  scales  are 
important  because  of  their  effect  on  the  large  scales.  Another  important  assumption  in  LES  is  that 
the  small  scales  of  turbulence  are  more  isotropic  and  more  universal  than  the  large  scales. 

The  usual  LES  starts  by  convoluting  the  Navier-Stokes  equations  (NSE)  with  a spatial  filter  gs- 

uH-  V • (uu)  - Re~*  Au  + Vp  = f , (1) 

where  u = gf  * u is  the  variable  of  interest.  The  filtered  NSE  (1)  do  not  form  a closed  system,  and 
most  LES  research  has  been  directed  at  modeling  the  subgrid-scale  (SGS)  stress 

T = uu  - u u.  (2) 

Arguably  the  most  popular  class  of  SGS  models  is  the  eddy  viscosity  type,  based  on  (variants 
of)  Sinagorinsky’s  model  [24], 

r = ~(Ca6)2  |Vsu|  Vsu,  (3) 
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where  <5  is  the  radius  of  the  filter  and  Vsu  = j(Vu  4-  Vu7 ).  Smagorinsky  proposed  (3)  with  C, 
a constant,  but  the  most  successful  models  have  been  those  in  which  Cs  is  computed  dynamically 
(the  dynamic  subgridscale  eddy  viscosity  model  [9]  and  its  variants). 

Other  successful  LES  models  include  the  scale-similarity  models,  the  mixed  models  (scale- 
similarity  models  coupled  with  an  eddy-viscosity  model  to  properly  dissipate  energy),  and  the 
RNG  models. 

Despite  its  undeniable  achievements,  LES  still  poses  numerous  challenges,  many  of  them  at  a 
very  fundamental  level.  One  of  these  challenges  is  the  need  for  more  mathematical  consistency,  which 
would  provide  the  means  for  sounder  mathematical  and  numerical  analysis  for  the  corresponding 
LES  models.  Other  fundamental  challenges  are  the  LES  boundary  conditions,  the  influence  of  the 
filter  on  the  SGS  model,  commutation  errors  in  filtering,  numerical  errors,  and  numerical  validation. 
A growing  current  of  opinion  in  the  LES  community  favors  the  need  for  new  approaches  that  could 
answer  these  fundamental  challenges. 

The  LES  model  used  in  this  paper  is  based  on  an  approximation  rather  than  a physical  analogy. 
This  model  has  evolved  in  several  steps.  First,  in  1974  Leonard  [16]  developed  a continuum  model 
of  u u 

S2 

u u = u u + ~A(u  u)  + OfiS4),  (4) 

where  S is  the  radius  of  the  filter  and  7 is  a parameter  in  the  Gaussian  filter.  Next,  in  1979  Clark, 
Ferziger,  and  Reynolds  [3]  developed  an  analogous  model  for  the  “cross  terms”  uu'  + u'u,  where 
u'  :=  u - u represents  the  turbulent  fluctuations.  The  approach  used  in  the  derivation  of  these 
models  involved  taking  the  Fourier  transform  of  the  corresponding  terms,  approximating  the  Fourier 
transform  of  the  Gaussian  filter  gg,  dropping  the  terms  of  order  0(S4)  and,  finally,  taking  the  inverse 
Fourier  transform  of  the  resulting  approximations.  Noticing  that  the  approximation  by  Taylor  series 
of  gs  actually  increases  the  high  wave  number  components,  Galdi  and  Layton  [8]  developed  a new 
LES  model  based  on  a rational  ((0,1)  Pade)  approximation  of  gs,  which  preserves  the  decay  of  the 
high  wave  number  components.  The  resulting  LES  model,  which  will  be  called  in  the  sequel  the 
“Rational”  LES  model,  reads  as  follows: 


w,  — Rc  1 Aw  4-  V - (ww)  4-  Vg  4-  V • 

w(x,  0)  = u0(x) 

Boundary  Conditions 


fS2 

— VwVw 

V2? 


) 


= f 


in  fl, 

in  D, 
on  90, 


where  (w ,q)  are  proposed  as  an  approximation  to  (u,p),  fi  C (d  = 2,3),  and 


tvwvw  ki-Esr 


1=1 


9w,  dwj 

dxi  9x; 


(5) 


A sound  mathematical  analysis  [2]  and  the  first  steps  for  a numerical  validation  [12]  were 
provided  for  the  “Rational”  LES  model  (5)  coupled  with  a Smagorinsky  model  with  a very  small 
coefficient  C,  (for  alternatives  to  the  Smagorinsky  model,  see  [14]). 
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This  paper  presents  numerical  results  for  the  “Rational”  LES  model  (5)  applied  to  the  3D 
channel  flow  test  problem  at  a Reynolds  number  based  on  the  wall  shear  velocity  Rer  = 180. 

2.  Numerical  Setting 

The  3D  channel  flow  (Figure  1)  is  one  of  the  most  popular  test  problems  for  the  investigation  of  the 
wall  bounded  turbulent  flows.  We  will  compare  our  numerical  results  with  the  extensive  data  exist- 
ing for  this  test  problem  (both  Direct  Numerical  Simulation  (DNS)  and  LES).  The  computational 
domain  is  periodic  in  the  streamwise  (x)  and  spanwise  (z)  directions,  and  the  pressure  gradient 
that  drives  the  flow  is  adjusted  dynamically  to  maintain  a constant  mass  flux  through  the  channel. 
We  selected  the  periodic  domain  sizes  (x  and  z)  as  Lx  — 27r,  Lz  — |tt,  and  we  chose  imity  as  the 
channel  half-width  {&chann'i  = !)•  We  employed  in  our  calculations  a 40  x 49  x 40  quadrature  grid. 
The  length  scale  8 is  computed  as  8 = (/A*  Az  A y(y),  where  Ax  and  A*  are  twice  the  largest 
spaces  between  the  Gauss-Lobatto-Legendre  (GLL)  points  in  the  x and  z directions,  respectively, 
and  Aj ,(y)  is  inhomogeneous  and  is  computed  as  an  interpolation  function  that  is  zero  at  the  wall 
and  is  twice  the  normal  meshsize  for  the  elements  near  the  wall  and  for  those  in  the  center  of  the 
channel.  For  the  “Rational"  LES  model  (5)  we  chose  7 = 6 (the  usual  choice  for  the  Gaussian 
filter). 

We  compared  the  “Rational”  LES  mode)  (5)  with  the  Smagorinsky  model  with  Van  Driest 
damping  (see  [18]).  The  Smagorinsky  model  with  Van  Driest  damping  is  obtained  by  modifying  the 
usual  Smagorinsky  model  (3)  to  reduce  the  eddy  viscosity  in  the  near-wall  region 

r = ~[Cs<f(l  - exp(-y+/A+))f  |Vsu|  V’u,  (6) 

where  y+  = (Schannci  — \y\)uT/v  is  the  nondimcnsional  distance  from  the  wall,  C,  is  the  Smagorinsky 
constant,  and  A+  is  a constant.  Here  uT  Ls  the  wall-shear  velocity,  uT  \fTwnu/p,  in  which  p is  the 
density,  and  tw„u  is  the  wall-shear.  In  our  calculations  we  made  the  usual  choices:  C,  = 0.1  and 
A+  = 25. 

The  numerical  simulations  wore  performed  using  a spectral  element  code  based  on  the  P/v  — P/v-2 
velocity  and  pressure  spaces  introduced  by  Maday  and  Patera  [17].  The  domain  is  decomposed  into 
8 x 10  x 8 elements,  as  shown  in  Figure  2.  The  velocity  is  continuous  across  element  interfaces  and 
is  represented  by  Nth-order  tensor-product  Lagrange  polynomials  based  on  the  GLL  points.  The 
pressure  is  discontinuous  and  is  represented  by  tensor-product,  polynomials  of  degree  N — 2.  Time 
stepping  is  based  on  an  operator-splitting  of  the  discrete  system,  which  leads  to  separate  convective, 
viscous,  and  pressure  subproblems  without,  the  need  for  ad-hoc  pressure  boundary  conditions.  A 
filter,  which  removes  5%  of  the  highest  velocity  mode,  is  used  to  stabilize  the  Galerkin  formulation 
[7].  Details  of  the  discretization  and  solution  algorithm  sue  given  in  [5],  [6]. 

Wc  applied  initial  conditions  consisting  of  the  parabolic  mean  flow  (Poiseuille  flow),  on  which 
a 2D  Tollmien-Schlichting  (TS)  mode  of  2%  amplitude  and  a 3D  TS  mode  of  1%  amplitude  were 
superimposed. 

In  our  numerical  experiments  we  considered,  as  a first,  step,  homogeneous  boundary  conditions 
for  the  “Rational”  LES  model  (5).  The  mathematical  and  numerical  investigation  of  more  realistic 
(slip-with-friction)  boundary  conditions  in  the  LES  context  began  in  [22],  [13]. 

The  numerical  results  include  plots  of  the  following  time-  and  plane-averaged  quantities:  the 
mean  streamwise  velocity  <K  u the  x,  y-component  of  the  Reynolds  stress  <S  uV  S>,  and  the 
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rras  values  of  the  streamwise  <SC  uV  wall-normal  v'v'  2>,  and  spanwise  w'w'  '3>  velocity 
fluctuations,  where  < • > denotes  plane  ( xz ) averaging,  <S  • S>  denotes  time  and  plane  (xz) 
averaging,  the  fluctuating  quantities  f are  calculated  as  f — /—  / >,  and  a “+”  superscript 

denotes  the  variable  in  waliunits. 

3.  Results  and  Conclusions 

These  numerical  results  were  the  basis  for  a posteriori  tests:  we  compared  them  with  the  fine  DNS 
simulation  of  Moser,  Kim  and  Mansour  [19]. 

Figure  3 shows  the  mean  streamwise  velocity;  note  the  almost  perfect  overlapping  of  the  results 
corresponding  to  the  two  methods  (“Rational”  LES  and  Smagorinsky  model  with  Van  Driest  damp- 
ing). The  same  behavior  can  be  observed  in  the  x,  y-component  of  the  Reynolds  stress  in  Figure  4, 
with  a slight  improvement  for  the  “Rational”  LES  model. 

Figures  5-7,  containing  the  rms  values  for  the  three  velocity  components,  merit  a more  detailed 
discussion. 

The  rms  values  of  the  streamwise  velocity  fluctuations  in  Figure  5 show  a better  (closer  to  the 
fine  DNS  benchmark  results  in  [19])  behavior  for  the  “Rational”  LES  model.  Near  the  wall  this 
improvement  is  clearer,  the  two  models  performing  similarly  away  from  the  wall. 

Figure  6 shows  the  rms  values  of  the  wall-normal  velocity  fluctuations.  Here,  the  Smagorinsky 
model  with  Van  Driest  damping  performs  slightly  better  near  the  wall,  and  consistently  better 
between  the  near-the-wall  region  and  the  center  of  the  channel.  Toward  the  center  of  the  channel, 
the  “Rational”  LES  model  performs  slightly  better. 

The  rms  values  of  the  spanwise  velocity  fluctuations  in  Figure  7 are  consistently  better  for  the 
“Rational”  LES  model,  with  the  exception  of  a very  short  portion  close  to  the  center  of  the  channel. 

In  conclusion,  the  “Rational"  LES  model  performs  better  near  the  wall  (with  the  exception  of 
the  rms  values  of  the  wall- normal  velocity  fluctuations),  whereas  the  Smagorinsky  model  with  Van 
Driest  damping  performs  better  closer  to  the  center  of  the  channel  (with  the  exception  of  the  rms 
values  of  the  spanwise  velocity  fluctuations).  The  two  models  are  comparable  in  accuracy. 

Further  research  should  be  directed  at  the  investigation  of  additional  terms  accounting  for  the 
discarded  0(SA)  terms  in  the  “Rational”  LES  model,  as  well  as  to  the  influence  of  the  length  scale 
S in  the  spectral  element  implementation  of  the  “Rational”  LES  model. 
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Figure  2.  Spectral  element  mesh. 


Figure  3 . Mean  streiunwise  velocity;  ...  fine  DNS,  o Sraagorinsky  with  Van  Driest  damping,  -f  “Rational”  LES 
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Figure  5.  rms  values  of  strcamwisc  velocity  fluctuations;  ...  fine  DNS,  o Smagorinsky  with  Van  Driest  damping,  -f 
“Rational*  LES. 
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Figure  6.  rm3  values  of  wall-normal  velocity  fluctuations;  ...  fine  DNS,  o Smagorinsky  with  Van  Driest  clamping,  + 
“Rational11  LES. 


Figure  7.  rms  values  of  spanwise  velocity  fluctuations;  ...  fine  DNS,  o Smagorinsky  with  Van  Driest  damping,  -f 
“Rational"  LES. 


